Structure of Neutron-Scattering Peak in both s ++ -wave and s±-wave states 

of an Iron pnictide Superconductor 
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We study the neutron scattering spectrum in iron pnictides based on the random-phase approxi- 
mation in the five-orbital model, for fully-gapped s-wave states with sign reversal (s±) and without 
sign reversal (s++). In the s++-wave state, we find that a prominent hump structure appears just 
above the spectral gap, by taking account of the quasiparticle damping 7 due to strong electron- 
electron correlation: As the superconductivity develops, the reduction in 7 gives rise to the large 
overshoot in the spectrum above the gap. The obtained hump structure looks similar to the reso- 
nance peak in the s±-wave state, although the height and weight of the peak in the latter state is 
much larger. In the present study, experimentally observed broad spectral peak in iron pnictides is 
naturally reproduced by assuming the s++-wave state. 

PACS numbers: 74.20.-z, 74.20.Rp, 78.70.Nx 



Since the discovery of superconductivity in iron pnic- 
tides with high transition temperature (T c ) next to high- 
Tc cuprates[l], the structure of the superconducting (SC) 
gap has been studied very intensively. The SC gap in 
many iron pnictides is fully-gapped and band-dependent, 
as shown by the penetration depth measurement [2] and 
the angle-resolved photoemission spectroscopy (ARPES) 
[3, 4], except for P-doped Bal22 [5]. The fully-gapped 
state is also supported by the rapid suppression in 1 jT\ 
(oc T n ; n ~ 4 - 6) below T c [6-8]. 

In iron pnictides, the nesting of the Fermi surface (FS) 
between hole- and electron-pockets is expected to in- 
duce the antiferromagnetic (AF) fluctuations in doped 
metal compounds. Since fully-gapped sign-reversing s- 
wave state (s±-wave state) is a natural candidate [9, 10], 
it is urgent to clarify the sign reversal in the SC gap via 
phase-sensitive experiments. One of the promising meth- 
ods is the neutron scattering measurement: Existence of 
the resonance peak at a nesting wavevector Q is a strong 
evidence for AF fluctuation mediated superconductors 
with sign reversal [11-13]. The resonance condition is 
lu ics < 2 A, where w rcs is the resonance energy and A is 
magnitude of the SC gap at T = 0. The resonance peak 
has been observed in many AF fluctuation mediated un- 
conventional superconductors, like high-Tc cuprates [14- 
16], CeCoIn 5 [17], and UPd 2 Al 3 [18]. 

Neutron scattering measurements for iron pnictides 
have been performed [19-22] after the theoretical pre- 
dictions [23, 24]. Although clear peak structure was ob- 
served in FeSeo.4Teo.6 [20] and BaFei.85Coo.isAs2 [21], 
its weight is much smaller than that in high-Tc cuprates 
and CeCoIns, and the resonance condition lu ics < 2 A is 
not surely confirmed, as we will discuss later. 

Nonmagnetic impurity effect also offers us useful 
phase-sensitive information. Theoretically, s±-wave state 
should be very fragile against impurities due to the inter- 
band scattering [25]; the predicted critical residual resis- 



tivity pf^,p for vanishing T c is about 20 /ificm. However, 
experimental pf^p reaches ~ 750 /ificm, which corre- 
sponds to the minimum metallic conductivity 4e 2 //i per 
layer [26]. Since this result supports a conventional s- 
wave state without sign reversal (s ++ -wave state), we 
have to resolve the discrepancy between neutron scatter- 
ing measurements and the impurity effects. 

In this letter, we study the dynamical spin suscepti- 
bility x s ( w iQ) based on the five-orbital model [9] for 
both S4.4. and s± wave states, and discuss by which pair- 
ing state the experimental results are reproducible. In 
the normal state, x s (u, Q) is strongly suppressed by the 
quasiparticle damping 7 due to strong correlation. How- 
ever, this suppression diminishes in the SC state since 7 
is reduced as the SC gap opens. For this reason, a promi- 
nent hump structure unrelated to the resonance mecha- 
nism appears in X s (w, Q) just above 2A in the s++-wave 
state. In the s±-wave state, very high and sharp reso- 
nance peak appears at u> les < 2A. We demonstrate that 
the broad spectral peak observed in iron pnictides is nat- 
urally reproduced based on the s++-wave state, rather 
than the s±-wave state. 

Now, we study the 10 x 10 Nambu BCS Hamiltonian 
Hk composed of the five-orbital tight-binding model and 
the band-diagonal SC gap introduced in rcf. [25]. The 
FSs are shown in Fig. 1 (a). Then, the 10 x 10 Green 
function is given by 



G(w„, k) 
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where Lo n = nT(2n + 1) is the fermion Matsubara fre- 
quency, G (F) is the 5x5 normal (anomalous) Green 
function, and is the self-energy in the ei-orbital basis. 
For a while, we assume that the SC gap for the a-th FS 
is band-independent; |A Q | = A. Hereafter, the unit of 
energy is eV, unless otherwise noted. 
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Here, we have to obtain the spin susceptibility as func- 
tion of real frequency. For this purpose, it is rather easy 
to use the Matsubara frequency method and the numer- 
ical analytic continuation (pade approximation). In the 
present study, however, we perform the analytical con- 
tinuation before numerical calculation in order to obtain 
more reliable results. The irreducible spin susceptibility 
in the singlet SC state is given by [13] 
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where x + = x + ui, k + = k + q, U = 1 ~ 5 represents the 
e?-orbital, and A (R) represents the advanced (retarded) 
Green function. pf v {x,k) = {G^,(x, k)-Gf v (x, k))/(2m) 

and p^\x,k) ee (F^ )A (x,k) - F^ )R (x,k))/(2m) are 



one particle spectral functions. Since p'f l ; F (x, k) = for 
|x| < A, Imx™,;,(o;,q) = for \u\ < 2A. That is, the 
particle-hole excitation gap is 2A. 

Then, the spin susceptibility x s ( w i<7) is given by the 
multiorbital random-phase-approximation (RPA) with 
the intraorbital Coulomb U, the interorbital Coulomb 
U', the Hund coupling J, and the pair- hopping J' [9]: 
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where vertex of spin channel Sf ia t t — U, U' , J and J' 
for l\ = ?2 = h = hi h = h 7^ h = h , h = h ^ h = h 
and l\ = Z4 7^ I2 = I3, respectively. Hereafter, we put 
J = J' = 0.15, U' = U - 2 J, U = 1 ~ 1.3, and fix the 
electron number as 6.1 (10% electron-doped case). In the 
present model, x s (0, q) takes the maximum value when q 
is the nesting vector Q — (7r,7r/16). Due to the nesting, 
X S (0,Q)/X°(0,Q) « 1/(1 - a S t) is enhanced; ast (< 1) 
is the maximum eig envalue of S°x m {0, Q) that is called 
the Stoner factor. 

In strongly correlated systems, x s ( uj tQ) is renormal- 
ized by the self-energy correction. In nearly AF met- 
als, for example, the temperature dependence of the self- 
energy induces the Curie- Weiss behavior of x s (0, Q)- At 
the moment, there is no experimental information on the 
k-, e-, and band-dependences of the self-energy. There- 
fore, we phenomenologically introduce a band-diagonal 
self-energy as Sj^(e) = ij(e)l- First, we estimate the 
value of 7(e) in the normal state. Since the conduc- 
tivity is given by a = e 2 J2 a N a (0)v^/2j(0), where 
N a (0) and v a are the density of states (DOS) and the 
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FIG. 1: (Color online) (a) FSs in iron pnictides. (b) in- 
dependence of Imx s (m,Q) for the s++-wave state (A = 0.4) 
and the normal state, where the unit of energy is eV. The 
"exact result" is obtained by eq. (2), and the "approximate 
result" is obtained by eq. (6). 



Fermi velocity of the a-th FS, we obtain p « (27[meV]) 
pttcm [25]. Since p(T) - p(0) ~ (5T[meV]) /xficm in 
BaFei.84Co .i 6 As2 below 100 K [27], 7(0) due to inelas- 
tic scattering is estimated as 2.5T which is comparable 
to that in over-doped cuprates. If we assume the relation 
7(e) oc (irT + e) in nearly AF Fermi liquid [28], we obtain 
7(e) ~2.5(T+e/7r). 

Now, we calculate Imx s (w,Q) in both normal and 
,s++-wave SC states, concentrating on the frequency 
uj ~ 2A. To estimate the renormalization of Imx s (u;, Q) 
due to the self-energy, we have to know the value of 7(e) 
with |e| ~ A in both normal and SC states. Considering 
that 7(e) = 2.5(T + c/tt) - 2A at T c = 2.2 meV and 
€ = A ~ 5 meV in BaFei.g5Coo.i5As2, in the present 
study, we simply put 7(e) in the normal state at T c as 



7(e) = 7o 



(4) 



with 70 > A. In the present model, ast = 0.84 (0.79) 
for U = 1.3 (1.2) when 7o = 0.1 and T = 0.002; the 
T-dependence of ast is small when 70 is fixed. 

In the SC state at T < T c , 7(e) = for |e| < 3A (= 
particle-hole excitation gap 2A plus one-particle gap A) 
[12], and its functional form is approximately the same 
as that in the normal state for |e| > 3 A. Then, we put 



7(e) = a(e)7 s 



(5) 



where (i) o(c) < 1 for |e| < 3A, (ii) a(e) = 1 for |e| > 4A, 
and (hi) linear extrapolation for 3 A < |e| < 4A. We have 
confirmed that the obtained results are insensitive to the 
boundary of |e (4 A in the present case) between (ii) and 
(hi). Although 7 S at T <C T c should be smaller than 70 
at T = T c , we simply put 7^ = 70 hereafter, which causes 
underestimation of the peak height of hxix s ■ 

Figure 1 shows Imx s (w,Q) obtained by eqs. (2) and 
(3) for U = 1.2, 70 = 0.4 and T = 0.01. In the s ++ - 
wave SC state, we put A = 70 and a(3A) = 0.05. In 
calculating eq. (2), we use 256 x 256 fc-meshes, and 
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1000 i-meshes. Although values of A and 7 in Fig. 1 
are very large to obtain enough numerical accuracy, the 
ratio 70/A ~ 1 is consistent with experiments. In the 
normal state, Imx s (u,Q) is suppressed by large quasi- 
particle damping 70 ~ A. In the SC state, the gap in 
lmx s (to,Q) is 2A. Since the particle or hole with en- 
ergy |e| < 3A is free from inelastic scattering in the SC 
state, the lifetime of particle-hole excitation with energy 
|e| < 4 A becomes long. For this reason, Imx s (w, q) shows 
a large hump structure for 2A < u) < 4A below T c in 
s ++ -wave state. 

Unfortunately, we cannot put smaller A and 7 in cal- 
culating eq. (2) in the five-orbital model, because of the 
computation time. To solve this problem, we perform the 
x-integration in eq. (2) approximately as follows: When 
7 = 7I, the retarded (advanced) 10 x 10 Green func- 
tion is expressed as Gm[m' ( x > *0 = E Q U^' a (x + (— )i-y - 
E%)- l U%' ' a * , where Eg (a = 1 ~ 10) is the eigenvalue of 
Hk and Uk is the corresponding 10 x 10 unitary matrix. 
We promise that = -E k +5 for 1 < a < 5. When 
7 is sufficiently small, then p l ,,(x,k) « ^ Q U k a S(x ~ 
' ' 5 - ),Q j and thus eq. (2) becomes 
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with IV.kq = 7 for 7 < 1. 

When 7 is as large as A, however, we have to check 
to what extent eq. (6) is reliable. Considering that the 
origin of the renormalization of x s is the quasiparticle 
damping i(E l k ) and j{E l k+q ), we introduce the following 
approximation: 



r„,, fe(} = 6-max{7(^) ) 7(4' +(7 )} 



(7) 



where b ~ 1 is a fitting parameter. f. q w in the 
SC state for \E l k \, \E l kJi _ \ < 3A, reflecting the absence 
of quasiparticle damping. In Fig. 1, we show numerical 
results given by the present approximation with b = 1.3; 
we replace 670 with 70 hereafter since b ~ 1. Since the 
"exact results" given by eq. (2) is quantitatively repro- 
duced, we decide to calculate Imx s (w,Q) using eqs. (6) 
and (7) for more realistic values of A and 7. We veri- 
fied that the present approximation works well when 7 is 
comparable to or smaller than A. 

Figure 2 shows Imx s (w,Q) obtained by eqs. (6) and 
(3) for U = 1.3 and T = 0.002. In the s ++ -wave SC state, 
we put A = 0.05; although it is a few times larger than 
the gap for Smllll with T c = 56K, it is enough smaller 
than the Fermi energies of electron- and hole-pockets 
[9]. In the numerical calculation, we use 1024 x 1024 k- 
meshes. Hereafter, we put a(3A) = O.OO3/70. When (a) 
70 = 0.003, Imx s (u> , Q) in the SC state approximately 
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FIG. 2: (Color online) Imx s (^,Q) for s++-wave (solid line) 
and normal (broken line) states, with 70 = ~ 0.1. 



equal to that in the normal state for to > 2A. As 70 
increases from (b) 0.05 to (d) 0.1, Imx s in the normal 
state decreases gradually, whereas that in the SC state 
depends on 70 only slightly, since 7(e) as for |e| < 3A. 
Therefore, in the case of 70 > A, Imx s (w, Q) in the SC 
state shows a prominent hump structure, and its peak 
value is about double of that in the normal state. In (d), 
experimental approximate "sum- rule" at fixed q = Q [21] 
is well satisfied. In Fig. 2 (c) and (d), a relatively large 
slope for |e| < 2 A is an artifact of the approximation due 
to large 70/A. 
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FIG. 3: (Color online) Imx s (i^,Q) for s++-wave (solid line) 
and normal (broken line) states, with A max = 0.07 and 
A min = 0.035. 

Next, we study the effect of band-dependent SC gap 
observed by ARPES measurements [3, 4]. In Fig. 3 
(a), we put U = 1.3, Ai. 2 , 4 = A max = 0.07eV for 
FS1,3,4, and A 2 = A min = 0.035eV for FS2. In (b), 
we introduce the anisotropy of the gap function for only 
FS3,4 with ratio 2; A& = A max (l — 0.5sin 2 #fc), where 
6 k = tan _1 (|fc y ( 2 .) 1/(1^(^)1 - 7r)) for FS3(4). Here, we 
put a(e) in eq. (5) as (i) O.OO3/70 for |e| < 3A min , 
(ii) 1 for |e| > 4A m i„, and (iii) linear extrapolation for 
3A min < |e| < 4A min . In Fig. 3 (a), Imx s (o;,Q) in- 
creases rapidly at ui = A max + A m j n = 0.105, and it 
shows a peak at w = 0.14. In (b), the peak is located 



at uj 



0.125, which is closer to A r 



0.105. 



In Fig. 3 (a) and (b), the width of the hump peak is 
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much sharper than that for the band-independent SC gap 
in Fig. 2, since luix s (w,Q) is reduced by damping for 
\ui\ > 4A min = 0.14. We have also calculated Iuix s (uj, Q) 
for A3 4 = A max and Ai .2 = A min , and verified that the 
obtained result is similar to Fig. 3. 

Here, we make comparison with experiments. The 
peak height and the weight in Fig. 3 (b) seems to 
be consistent with the neutron scattering measurements 
in iron pnictides [19-22]. In BaFei.gsCoo.isAss (T c — 
25K), the observed "resonance energy" is u; ros = 9.5 
meV [21]. According to ref. [3], A max /T c w 3.5 and 
A m i n /A max w 0.35 in many iron pnictides. (More smaller 
A max , min is reported in ref. [2].) Thus, A max + A min « 
4.7T C = 10 meV is comparable to w ros = 9.5 meV 
in BaFei.85Coo.i5As2. Moreover, finite Imx s (u,Q) for 
u) > 0.3u> Tes in ref. [21] may suggest the existence of 
SC gap anisotropy. Therefore, the theoretical result in 
Fig. 3 (b) is well consistent with experimental data. We 
have verified that the hump structure of Imx s (w, q) with 
q = (tt, 0) is very small for 70 ~ 0.1. 
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fermion Kondo insulator CeNiSn. As shown in Fig. 1 of 
Ref. [31], neutron scattering spectrum at q = (0,7r,0) in 
CeNiSn shows a prominent hump peak structure above 
the hybridization gap below the Kondo temperature 7k, 
which looks very similar to the spectrum observed in 
iron pnictides below T c [19-22]. This hump structure 
is well reproduced by the dynamical-mean-field theory 
based on the periodic Anderson model [32]. This fact 
demonstrates that large hump in lmx s (ui,Q) can ap- 
pear in strongly correlated systems with one-particle gap, 
without the necessity of the resonance mechanism. 

In summary, we have studied Imx s (w, Q) in iron pnic- 
tides based on the five-orbital model, and revealed that 
a prominent hump structure appears just above 2A in 
the s++-wave state, by taking the strongly correlation 
effect via 7. This hump structure becomes small as a s 
decreases in the over-doped region or q deviates from the 
nesting. At present, experimental data can be explained 
in terms of the s++-wave state very well. Further ex- 
perimental efforts are required to determine the height 
and width of the "resonance peak", and the magnitude 
relation between uj rcs and A max + A m ; n . 

This study has been supported by Grants-in-Aid for 
Scientific Research from MEXT of Japan, and by JST, 
TRIP. Numerical calculations were performed using the 
facilities of the supercomputer center, ISSP. 
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FIG. 4: (Color online) Imx s (u,Q) for s±-wave (solid line) 
and normal (broken line) states, with U = 1.2 and 1.0. 

We also analyze the resonance peak for the s-j--wave 
state in Fig. 4. In this case, the spin wave without 
damping is observed as the "resonance peak" at uj ics < 
2 A. Figure 4 shows the numerical results for (a) U = 1.2 
and (b) U = 1.0 in the case of A1.2 = -A 3 , 4 = 0.05. 
In (a), a very sharp and high resonance peak appears at 
<^rcs = 0.85 < 2A, consistently with previous theoretical 
studies [23, 24]. The case (b) with U = 1.0 corresponds 
to the "heavily overdoped" since ast = 0.69 and T c ~ 0. 
The obtained resonance peak in Fig. 4 by taking 7(e) 
into account is too large to explain experiments even in 
the case of agt = 0.69. In Bi-based high-T c cuprates, the 
width of the resonance peak is wide due to the sample 
inhomogeneity (i.e., nanoscale distribution of T c ) [16]. 
However, weight of the peak is 10 times larger than that 
in BaFei.85Co .i 5 As2 [21]. 

In the present study, we have neglected the impurity 
effect since its influence on x s { UJ ^Q) is expected to be 
small. In fact, in the single band model, the reduction in 
X° due to the impurity self-energy is almost canceled by 
the impurity vertex correction [29]. Moreover, impurity 
effect tends to enhance x s { ljJ iQ) m the modified FLEX 
approximation in nearly AF metals [30] . 

Before closing the study, we shortly discuss the heavy 
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